Does Silicene on Ag(lll) Have a Dirac Cone? 
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We investigate the currently debated issue of the existence of the Dirac cone in silicene on an 
Ag(lll) surface, using first-principles calculations based on density functional theory to obtain the 
band structure. By unfolding the band structure in the Brillouin zone of a supercell to that of 
a primitive cell, followed by projecting onto Ag and silicene subsystems, we demonstrate that the 
Dirac cone in silicene on Ag(lll) is destroyed. Our results clearly indicate that the linear dispersions 
observed in both angular-resolved photoemission spectroscopy (ARPES) [P. Vogt et al, Phys. Rev. 
Lett. 108, 155501 (2012)] and scanning tunneling spectroscopy (STS) [L. Chen et at, Phys. Rev. 
Lett. 109, 056804 (2012)] come from the Ag substrate and not from silicene. 



Silicene, a two-dimensional (2D) allotrope of Si with a 
hexagonal honeycomb-like lattice similar to graphene [l[ , 
has recently attracted intense attention. It was proven in 
first-principles studies that low-buckled silicene is ther- 
mally stable and has a linear electronic dispersion near K 
points at corners of the first Brillouin zone 0411 , which is 
similar to the behavior of graphene. As a result, Si atoms 
in silicene are conjectured to be partially sp^ hybridized 
0. Silicon nanostructures and silicone were successfully 
synthesized by depositing Si atoms on surfaces of Ag 0- 
li], ZrBa and recently Ir 0. The Ag(lll) sub- 
strate is ideal for growing silicene because the tendency 
to form an Ag-Si alloy is low Q, and as a result there 
are several different atomic arrangements for silicene on 
Ag(lll) surfaces. A linear dispersion relation was found 



in the 4x4 structure by Fleurence et al. [15|. Feng 
et al. reported three different phases of silicene on 
Ag(lll) surface, two with a 4 x 4 unit cell with respect to 
the Ag(lll) lattice, and a third phase, a \/3 x -s/S recon- 
struction in reference to the low buckled silicene lattice, 
in which they also found the existence of Dirac fermions 
0. Atomic arrangements a/13 x -\/T3 R 13.9° of silicene 
on an Ag(lll) surface were also observed [l^l- Various 
configurations were constructed and simulated numeri- 
cally using density functional theory (DFT) 
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A common motivation of these studies was to detect 
and utilize the Dirac fermions in silicene. Two groups 
have claimed to find evidence for the existence of Dirac 
fermions in silicene on Ag(lll) surface. In one exper- 
iment 0, a linear dispersion near the Fermi energy in 
the silicene- Ag( 111) system was observed using angular- 
resolved photoemission spectroscopy (ARPES). No such 
band was observed for a bare Ag(lll) surface, and so the 
linear dispersion was attributed to silicene. The Dirac 
point was measured to be 0.3 eV below the Fermi energy, 
and the Fermi velocity was estimated to be 1.3x 10^ ms~^ 
0. In a second experiment the quasiparticle inter- 
ference (QPI) patterns at the surface of the silicene- 
Ag(lll) system were observed using scanning tunnel- 
ing spectroscopy (STS). The Dirac point was deduced 
from the linear dispersion curve to be ~ 0.5 eV below the 
Fermi energy, and the Fermi velocity is ^ 1.2 x lO^ms"^. 



The Fermi velocities from these two experiments are 
close to the theoretical prediction of ^ 10^ ms^^ Q. If 
judged merely from Fermi velocities, the experimentally 
observed linear dispersions do coincide with the theoret- 
ical band structure. However, the energy windows where 
the two experiments 0, 01 observed a linear dispersion 
arc too large. Band structure from first-principle calcu- 
lations showed linear dispersion only within a energy 
interval ±0.4 eV, which is several times smaller than the 
reported energy ranges for linear dispersions, — 3.0eV to 
-0.3 eV and 0.4 eV to 1.2 eV 0. This discrepancy 
and the lack of dispersion measurements in the vicinity 
of and across the Dirac point make the evidence for Dirac 
fermions in silicene rep orted by Ref. inconclusive. A 
very recent paper [l7| reported Landau level measure- 
ments in silicene on Ag together with band structure cal- 
culations. Absence of characteristic signals attributed 
to the Landau levels disagrees with the experimental re- 
ports @, 01 ■ Band structure calculations showed that 
Pz orbitals are strongly hybridized with and delocalized 
into the Ag substrate. The authors argued that the lin- 
ear dispersion observed in Ref. @ is not contributed by 
silicene, and the only reasonable explanation is that it 
comes from the Ag substrate 17|; however, there is no 



clear and straight-forward evidence presented. 

In this Letter, we report results from DFT calcula- 
tions that aim to understand the electronic structure of 
silicene on an Ag surface in the context of the linear dis- 
persion observed in experiments @, 0| • In all the above- 
mentioned work, DFT was routinely used to construct 
an atomic structural model that reproduces the observed 
STM images [1, 0, [H - 15 1 . Because of distortions of sil- 
icene on Ag surface, the unit cell of the silicene- Ag calcu- 
lation is 3 X 3 times large as primitive cell of silicene. The 
resulting band foldings make the calculated band struc- 
tures in Ref. 17[ too complex to extract useful informa- 
tion. Although Ref. [l3l argued that the experimentally 
observed linear dispersions are not from silicene but from 
the Ag surface, no direct connection was made between 
the experimental linear dispersions and Ag bands. In 
this work, we introduce a modified effective band struc- 
ture (EBS) technique to unfold bands from supercell cal- 
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culations and thus make tractable the identifying of sig- 
nals from experiments with specific bands. Unfolding the 
bands enables us to clearly identify the origins of the lin- 
ear dispersions reported in Refs. Q and Q. 

We used the Vienna Ab initio Simulation Package 
(VASP) [l^, a code based on a plane- wave basis and 
projector-augmented wave datascts [l8|. The Perdew- 
Burker-Ernzerhof form [l^ of exchange-correlation func- 
tional was used in all the calculations. The calculated 
bulk lattice constant oi fee Ag is 4.174 A. Although the 
atomic arrangements of silicene on Ag(lll) surfaces are 



not unique [a, [12, Il3| , we focused on the 4x4 atomic 



arrangement since it is the most common structure [l4l |. 
The geometric configuration of silicene on an Ag(lll) 
surface with a 4 x 4 atomic arrangement was constructed 
and optimized (see Fig. [2a)(b)). In one unit cell of the 
slab used for surface studies there arc five Ag layers, each 
with (4 X 4) Ag atoms. Ag atoms in the bottom two 
layers were kept fixed at bulk lattice positions. On the 
Ag(lll) surface, the silicene sub-unit consists of 18 Si 
atoms. After relaxation (with forces on atoms smaller 
than 0.01 eV/A), positions of six of the 18 Si atoms shift 
upward (away from Ag) with respect to the other twelve, 
and the average distance between the silicene sheet and 
the Ag surface is 2.44 A, which is close to that reported 
in Ref. Geometric relaxation was also performed for 
standalone, low-buckled silicene. 

For optimized geometries we then calculated the den- 
sity of states (DOS). For standalone low-buckled silicene, 
the DOS at the Fermi energy is zero (see Fig. HJc)). 
Our symmetry analysis of wave functions showed that in- 
deed at the Fermi energy the 7r-bonding and anti-bonding 
bands do not overlap but touch each other at K points at 
corners of the first Brillouin zone, indicating that stan- 
dalone silicene is a zero-gap semiconductor. When sil- 
icene is placed on the Ag surface, the DOS projected 
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FIG. 1. (color online) (a) Top and (b) side view of the 
silicene- Ag system (small balls: Si; big balls: Ag). (c) (pro- 
jected) density of states of low buckled silicene (line) and of 
silicene on Ag(lll) surface (circles). 



onto the silicene is no longer zero at the Fermi energy (see 
Fig.[ljc)); instead, a peak emerges near the Fermi energy 
as a result of silicene- Ag interaction. The net charges on 
Ag and Si atoms are negligible (less than 0.03 e~ per 
atom according to Bader charge analysis [l^l), but the 
difference between the DOS of pristine silicene and the 
DOS projected on silicene can not be describe by a sim- 
ple rigid shift; the interaction between Ag substrate and 
silicene is therefore beyond charge doping. 

Next, we looked into band structure, which offers more 
information on electronic structure than the DOS, to 
compare the band structure of the silicene-Ag(lll) sys- 
tem with that of standalone low-buckled silicene. How- 
ever, due to the distortion of silicene when deposited on 
Ag surfaces @, the unit cell is 3 x 3 times the primitive 
unit cell of low- buckled silicene Q . As a result of band 
foldings, a direct comparison between the band structure 
from supercell calculations and that from primitive unit 
cell calculations is meaningless [l3|- Even if one per- 
forms a supercell calculation for the standalone silicene, 
it is not trivial (nearly impossible) to extract informa- 
tion that can be compared with experimental data; band 
unfolding is necessary. 

There are several options to unfold the band structure 
from supercell calculations. First is to use the orbital- 
resolved spectral functions that can be calculated using 



the Wannier function based unfolding method l2l| . This 



method is well suited to our system, but the Wannier 
functions used in this method make it complicated to 
apply. A second possibility that is free from the compli- 
cations involving Wannier functions is to use the effective 
band structure (EBS) method 2^ 2^- This method was 
originally used to study the effective dispersion in alloys; 
however its original form is not ready to be used for our 
system because of its lack of atom and orbital resolutions. 
We thus introduced a modified EBS method, which en- 
ables calculations of orbital-projected spectral functions. 

In the original EBS method [2^, the spectral function 
is defined as 



(1) 



JV 



where the j^{k) are inner products of Kohn-Sham 
(KS) wave functions \KjN^ in the first Brillouin zone 
(BZ) of the supercell and KS wave functions 1^:^71^ in the 
first BZ of the primitive cell; E^^^ arc cigcnencrgies on 

predefined fc-points {Kj} in the first BZ of the supercell 
system. 

To extend Eq. ([1]), we define the projected spectral 
functions on specified atomic orbitals |a) as 
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(a) Si-projected EBS 



(b) Ag-projected EBS 




FIG. 2. (Color online) Si- and Ag-projected effective band structure (EBS) of the silicene-Ag system, (a) Si-projected 
effective band structures, with band structure of standalone low-buckled silicene plotted as blue dashed lines. Note that 
the band structure of low-buckled silicene was shifted downward by 1.1 eV. (b) Ag-projected effective band structures with 
experimental observed linear dispersions (red circles: measurements from Ref. [y|, blue squares: measurements from Ref. [3|) 



where the Pa^'^ are projection functions defined as 

p§^^ = {K,N\a) 



(3) 



i.e., expansion coefficients of KS wave functions P^^ j^{k) 

in atomic orbitals |q;). With this extension, we calcu- 
lated the spectral functions projected onto Ag and Si 
orbitals with k along a high-symmetry path T-K-M-T in 
the first BZ corresponding to the unit cell of standalone 
low-buckled silicene. The calculated spectral functions 
were then summed into Ag- and Si-projected spectral 
functions and divided by the number of atoms of a kind 
in one unit cell, so their magnitudes are ready to com- 
pare. These spectral functions are functions of k and E; 
magnitudes of the spectral functions at (fc, E) are rep- 
resented in Fig. [3] by different gray levels (darker points 
represent larger spectral magnitude at that point). Peaks 
in these spectral functions constitute continuous bands. 

In Fig. IHa), the band structure of standalone low- 
buckled silicene is also plotted with the Si-projected spec- 
tral function as the blue dashed lines. The band structure 
of standalone silicene was shifted 1.1 eV downward in or- 
der to match the Si-projected spectral functions around 
— l.leV, which suggests the interaction between Ag sur- 
face and silicene is more than charge doping. The Ag 
d-bands dominate in the energy interval — 6eV through 
— 3eV, so the Si-projected spectral function are hardly 
visible in this energy range (see Fig. Elja)). Most bands 
appear in both the Si- and Ag-projected spectral func- 
tions, because the corresponding wave functions arc nei- 



ther entirely localized in the Ag substrate nor in the 
silicene sheet. Magnitudes of the Si- and Ag-projected 
spectral functions can be compared in order to tell which 
atoms these bands come from. The projected EBS shown 
in Fig. [2] are somewhat complex, because we use in our 
calculations a slab with finite thickness, and the contri- 
butions from surface states mix with bulk states [24 1. 

If one plots the band structure for a system with a 
Dirac cone along high-symmetry path, the Dirac cone has 
two branches touching each other at the K point (Dirac 
point). It seems from Fig. [2][a) that the Dirac point in 
silicene on the Ag substrate moves to — l.leV, because 
the upper branch of the Dirac cone in the T-K portion 
and the lower branch in the K-M portion coincide with 
that of standalone low-buckled silicene moved downward 
by — l.leV. However, compared to Fig. [5Jb), the lower 
branch of the Dirac cone in the T-K portion is not dom- 
inated by Si. Most important, the upper branch of the 
Dirac cone in the K-M portion no longer exists. These 
observations indicates that the Dirac cone in silicene on 
Ag has been destroyed. 

The distortion of silicene on an Ag(lll) surface and 
its interaction with Ag(lll) surface are two possible rea- 
sons for the disappearance of Dirac cone. In order to 
show how the distortion affects the band structure of sil- 
icene, we did another EBS calculation in which all the Ag 
atoms in the unit cell were deleted, while silicon atoms 
were kept at their positions. The effective band struc- 
ture and density of states are shown in Fig. |31 The most 
import observation from Fig. |3] is the ^ 0.3 cV energy 
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FIG. 3. (Color online) (a) Effective band structures and 
(b) density of states of standalone silicene with additional 
buckling are compared with that of low-buckled silicene (red 
dashed lines). 



persion observed in Ref. [6| is —3 eV. Overall, based on 
our calculations we reach a different understanding of the 
experimental observed linear dispersions: they are from 
Ag(lll) surface instead of from silicene. 

In summary, inspired by experimental observations of 
linear dispersions in silicene-Ag(lll), we conducted band 
structure calculations using DFT and calculated the pro- 
jected effective band structures. The Dirac cone in sil- 
icene was observed to be destroyed by distortions and 
interactions with Ag surface. The linear dispersions ob- 
served in experiments were found to coincide with the 
Ag-projected but not Si-projected EBS, which indicates 
that these linear dispersions come from the Ag surface in- 
stead of from silicene. We presented a different interpre- 
tation of the experimental linear dispersions, attributed 
to be from the Ag surface instead of from silicene. 

This work was supported by the US Department of 
Energy (DOE), Office of Basic Energy Sciences (BES), 
under Contract No. DE-FG02-02ER45995. The compu- 
tation was done using the utilities of the National En- 
ergy Research Scientific Computing Center (NERSC). 
Authors also would like to thank T. Berlijn for fruitful 
discussion. 



gap opening at the K point. The 7r-bonding and anti- 
bonding bands near the Fermi energy do not touch each 
other any more, and develop to two flat bands in the K-M 
path which correspond to the two peaks near the Fermi 
energy in the DOS plot. Our conclusion from this calcu- 
lation is that the distortion seen for silicene on Ag(lll) 
surface alone can destroy the Dirac cone in silicone. 

Finally, we discuss the comparison between our results 
and experimental observations. Based on the (projected) 
EBS, we have concluded that the Dirac cone in silicene 
is destroyed by the distortion of silicene and the inter- 
action with Ag surface. However, linear dispersions over 
very large energy ranges were observed in experiments 
d, 0]. To understand this, we plotted the linear dis- 
persions observed in both experiments together with our 
Si-projected EBS of the silicene- Ag( 111) system, but we 
found no coincidence between the theory and the experi- 
ments. We do find a connection when placing the exper- 
imental dispersion functions on the Ag-projected EBS 
(see Fig. [IJb)): we observe that the experimental dis- 
persions coincide with the sp-band of the Ag(lll) sur- 
face. Although sp-electrons in Ag are free-electron-like 
and the sp-bands of Ag are actually parabolic, fitting to 
a linear function is acceptable for a finite energy interval 
away from its minimum energy. The large energy range 
where the sp-band of Ag appears is consistent with that 
in which the experiments observed linear dispersions, and 
provides the only plausible explanation. Additionally, 
below —3 eV, the sp-bands of Ag are buried by d-bands, 
which can explain why the lower bound of the linear dis- 
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